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ABSTRACT 

With the observations of the CS 5-4 line toward a sample of 24 infrared bright 
galaxies using HHSMT, we detected CS 5-4 emission in 14 galaxies, including 12 
ULIRGs/LIRGs and 2 nearby normal galaxies. As a good dense gas tracer, which has 
been well used for studying star formation in the Milky Way, CS 5-4 can trace the 
active star forming gas in galaxies. The correlation between CS 5-4 luminosity, which 
is estimated with detected CS 5-4 line emission, and infrared luminosity in these 14 
galaxies is fitted with correlation coefficient of 0.939 and the slope close to unity. This 
correlation confirmed that dense gas, which is closely linked to star formation, is very 
important for understanding star formation in galaxies. 
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1 INTRODUCTION 

The star formation law, the relation between star formation 
rate (SFR) and local gas density, is a fundamental tool to 
study star formation processes and underst and galaxy evo - 
lution. This relation was first proposed by ISchmidtJ (jl959n 



as psfr ~ (pgas) , where psfr and p gas are the vol- 
ume densities of SFR and gas, respectively, while N ~ 2 
and the atomic gas density was used for discussion. With 
measurements of disk-averaged total gas surface density 
E gas = Ehi + Eh 2 and the disk-averaged star formation 
surface density Esfr. in 61 normal spiral galaxies a nd 36 
infrared-selected starburst galaxies, iKennicuttl (1 19981 ) gave 
the best fit for the Schmidt law with index N — 1.4 ± 0.15. 
On the other hand, with pixel-by-pixel analysis at sub- 
kpc scales of 18 nearby galaxies, whi ch include 7 spira l 
galaxies and 11 late- type/dwarf galaxies. iBigiel et al.l (|200cf ) 
found that a best-fit molecular Schmidt law of Esfr = 



10 



-2.1±0.2^1.0±0.2 



Similar results (N ~ 1) were ob tained 



with the data of 23 nearby galaxies (|Lerov et al.ll2008T ). while 
the slope changes to 1.17 usi ng galaxies includin g z ~ 1 — 3 
normal star- forming galaxies (jGenzel et alJlioiol ). But if the 
ultraluminous infrared galaxies (ULIRGs) with extremely 
high star formation rate, which have the star formation ef- 
ficiency four to 10 times t han those normal star forming 
galaxies (jGenzel et alj|2010h . are included, the slope will be 
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higher tha n the results fr om Bigiel et al.l (|2008l ) , iLerov et al.l 
i|200Sf ) and lGenzel et al.l <|2010h . 

Lines of molecules with high dipole momentum, such as 
HCN, HNC, CN, HCO+ and CS, are thought to be better 
tracers of dense star forming gas than CO lines. These lines 
have been used to study star formation in the Milky Way, 
while they are generally difficult to be detected in extra- 
galactic object s due to their faint emission . HCN, HNC and 
HCO + lines l|Gracia-Carpio et al.l 120081 : iBussmann et al.l 
have been detected in a few galaxies with millime- 
ter (IRAM 30m) and sub-millimeter t elescopes (JCMT , 
HHSMT). CN line in several galaxies (|Aalto et al.1 12002T ) 
also has been observed with SEST. With the observations 
of den se gas tracers toward 65 galaxies (|Gao fc Solomon! 
l2004ah . which is the largest extra- galatic sample wi t h dense 
gas tracer observation up to now, Gap fc Solomon! (|2004bl ) 
found a strong linear correlation between HCN and IR lumi- 
nosity (N ~1) with more than 3 orders of infrared luminosity 
range, including LIRGs and ULIRGs. T his relation has been 
extended to the Milky Way dense cores dWu et al.ll2005l) and 
possibly high-z galaxies and QSOs as well ( Gap et akl 20071 ) . 
Such results agree well with that dense melecular gas is more 
directly related to final star formation than total molecular 
gas traced by CO lines from local universe to high-z. 

The line ratios of HCN, HNC, CN, and HCO+ molecules 
vary in different galaxies due to different chemical and phys- 
ical conditions, and also maybe due to pumping mechanism, 
which means the dense gas mass estimated with the obser- 
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vation of only one dense gas tracer may be biased. As a good 
dense gas tracer, CS has been widely used to study dense 
cores in GMCs in the Milky Way. Since it is not affected 
by shock chemistry, which is the case for HCO + molecule, 
and it is not related to HCN, HNC, and CN chemistry, 
CS can give an independed view of dense gas in galaxies 
from those based on HCN, HNC, CN, and HCO+ With 
critic al density of about 10 times that of HCN 1-0 (|Evansl 
Il999h . CS 5-4 can trace dense gas in active star forming re- 
gions. CS 5-4 line has been d etected in three nearby galaxie s 
(IC342, NGC 253, and M82) (|Mauersberger fe Henkejl989h . 
one nea rby major merger NGC 4038/4039 (the Antennae 
galaxy) (iBavet et alj|200l ). and the nearest ULIRG Arp 220 
(|Greve et alj|2009h with JCMT. CS 5-4 in several galaxies, 



including sources have been observed with other telescopes 
listed before and NGC 1068 fc Maffei 2 have been detected 
with IRAM 30M (|Martm et al.1 120091 '). Most of such galax- 
ies are not ULIRGs/LIRGs, which are galaxies with the 
most active star formation in local universe. So with ob- 
servations of CS 5-4 toward a sample of galaxies with more 
total sources and including more ULIRGs/LIRGs, we can 
see if dense gas traced by different molecules can give similar 
correlation or not, by co mparing with existing HCN results 
(Gao fc Solomon! l2004bh . 

In this letter, we describe the observations and data 
reduction in §2, present our main results and discussion in 
§3, and give a brief summary in §4. 



2 OBSERVATIONS AND DATA REDUCTION 

Our sample is composed of 24 galaxies including 20 local 
(U)LIRGs and 4 nearby normal galaxies that have strong 
CO emissions. The (U)LIRGs are selected from RBGS 
(jSanders et al.l l2003l l with S > -20°, L IR > 10 11 L , 
Fioofim > 25 Jy, Fao^m > 20 Jy. The nearby galaxies in our 
sample ha ve been observed with millimeter telescopes for 
HCN 1-0 (|Gao fc Solomonll2004bh . which provides enough 
dense molecular gas traced by HCN 1-0, for higher possi- 
bility of CS 5-4 detection. Table 1 listed the information of 
sources in this sample, including source name, the coordi- 
nates (J2000), luminosity distance, red-shift, and also the 
infrared luminosity. The sources in Table 1 are ordered with 
the luminosity distance. 



2.1 CS 5-4 observations with HHSMT 

The 10 m Heinrich Hertz Submillimeter Telescope 
(HHSMT) on Mt. Graham, Arizona, was used to observe the 
CS 5-4 line in this sample. These observations were done in 
good (t ^ 0.3) to medium (r ~ 0.3-0.6) weather conditions, 
with typical system temperature of 210 K, during several 
runs on the HHSMT in February 2009, December 2009, and 
February 2010. The 1.3 mm ALMA receiver was employed 
with lower side-band tuning to the CS 5-4 frequency in the 
single-sideband (SSB) mode. The beam size at the frequency 
for CS 5-4 (~ 2UGHz) is ~ 33". In order to ensure the 
baselines as flat as possible, most observations were done in 
fast beam-switching mode with a chop throw of 2' in az- 
imuth (AZ) in a chopping frequency of 2.2 Hz. For a few 
nearby galaxies, position switching mode was used to avoid 
the background contamination. During the 2010 February 



run, the ALMA receiver had only H-pol working properly, 
while the dual polarization mode was used on other time. 

We used both the Acoustic Optical Spectrometer (AOS) 
and Forbes Filter Banks (FFB) simultaneously, with ~ 1 
GHz bandwidth and 1024 channels to cover the full velocity- 
width of the expected broad lines. The data of both spec- 
trometers were used to check with each other, but not com- 
bined together. Most of the final spectra were used the FFB 
data, for its relative stability. Our observing strategy in- 
volved pointing and calibration observations every 2 hr using 
Saturn, Jupiter or Venus when available. Calibration scans 
were obtained in position-switching mode with reference po- 
sition of 5' in RA. We found typical pointing errors of 2-3" 
and measured the main-beam efficiency, rj m b, to be 0.6 for 
filter-bank A (H-polarization) and 0.8 for filter-bank B(V- 
polarization). The systematic flux calibration uncertainty is 
about 20%. 



2.2 CS data reduction 

All CS data reduction was performed using the CLASS pro- 
gram of the GILD AS package developed at IRAM. The ve- 
locity intervals based on the HCN or CO line profiles in 
literature are used to set the velocity range for all sources. 
Each baseline was inspected by eye and quantified into dif- 
ferent levels by their baseline flatness, system temperature, 
standing wave, and etc. About 10% -15% bad spectra were 
discarded from the qualification. For each source, the spectra 
were co-added, weighted by the rms noise of each spectrum. 
The final spectra are zoomed in the velocity range to a lit- 
tle bit wider than CO 1-0 line width from literature, and a 
zero level baseline is subtracted to get the 'local' baseline 
level. We smoothed the spectra to resolution from ~20 to 
50 km s _1 , depending on the noise levels, and convert from 
an antenna temperature, Ta, to the main-beam tempera- 
ture, T m t, by scaling the main-beam efficiency, rjmb, using 
T-mb = T a* /rj m b- Then the CS 5-4 line flux in each galaxy 
was integrated with the velocity range of CO line from liter- 
ature, and listed as the second column in Table 2, including 
the integrated flux and error bar. 

Single component gaussian profile is also used to fit each 
of the spectra with solid detections, to derive the area, sys- 
tem velocity, and line width of these lines. The fitted line 
flux of each galaxy was listed in the third column, while the 
fourth column was the line width. 

The upper limits of the CS measurements correspond to 
sources without solid detections, where 3 a are used in the 
following analysis. The upper limits were calculated through 
the relation 

orprms A 

T a Av (Kkms- 1 ) a (1) 
sj Av/ Av rcs 

where Av is the FWZI (Full Width Zero Intensity) line 
width of galaxies; Av reB is the velocity resolution (in our 
case « 20 — 50 kms -1 after smoothing) and T™ s is from 
the baseline fitting. 



2.3 CS line luminosity and the infrared 
luminosity correction 

Most of the dense gas emission in galaxies with distance 
larger than that of NGC 2276 would be expected to concen- 
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trate within the SMT beam (~ 33" for CS 5-4). We would 
like to consider that we have picked up most of the CS 5-4 
emission and will compare with the total infrared luminos- 
ity of entire galaxy in those sources. On the other hand, 
infrared luminosity in 6 sources from NGC 6946 to NGC 
2276 at the beginning of Table 1 have been corrected from 
the whole galaxy to the region within SMT beam. We used 
MIPS 24/^m image of those galaxies to obtain the ratio of 
such region to the whole galaxy and then derived infrared 
luminosity with CS 5-4 beam b y scaling the total in frared 
luminosity listed in Table 1 from ISanders et all ((2003) with 
24/^m flux ratio. 

We computed the CS line luminosity using 

L' CS = f T b {CS)dVdQ, a d 2 A ^ ■K/(4ln2)ef nb I CS dl(l+z)- 3 (2) 

, where Ics is defined as: 

Ics= fT mb {GS)dV= f T b (CS)dVdn s /[(l + z)n s * b ] (3) 

, where L' cs is in the unit of K km s~ 1 pc 2 . These 2 f ormula e 
were the formulae (1) and (2) in I Gap fc Solomon! (|2004al l . 
where we used CS 5-4 instead of HCN 1-0. All of the mea- 
sured line fluxes along with their errors are tabulated in 
Table 2. The errors include thermal and calibration errors, 
and also errors on the assumed aperture efficiencies. Each of 
these is assumed to be the uncorrelated and thus is added 
in the quadrature. The noise level was calculated for each 
line using RMSu ne = RMSchannei x VN channe i x FWZI, 
assuming all channels have the same response and are with- 
out any intereference between each other. 



3 RESULTS AND DISCUSSION 

We detected CS 5-4 with velocity integrated flux higher than 
3(J level in 13 galaxies. We used the CO line width to obtain 
velocity integrated flux, exc ept for NGC 6240 i n which we 
used CS 7-6 line range from iGreve et al.l (|2009h . Since CO 
line width is wider than that of dense gas tracers such as 
CS 5-4, we may over-estimate the noise level and all the 13 
sources with velocity integrated flux higher than 3a should 
be solid detections. We show spectra of 2 detected galax- 
ies in the top (Arp 220) and bottom (NGC 1068) of Fig- 
ure 1. CS 2- 1, 3-2, 5-4 and 7-6 lines in Arp 220 have been 
observed by IGreve et al.l (|2009l ) and showed different line 
shapes for different transitions. The velocity range and to- 
tal intensit y of our CS 5-4 sp ectrum agree well with the 
result from IGreve et all (|2009h . but with a little bit differ- 
ence in the line shape. Our CS 5-4 spectrum in NGC 6240 
sh ows a similar veloc ity range to that of CS 7-6 observed 
by IGreve et all (|2009h . which is about 200-300 km/s blue 
shifted than CO lines. Even though the integrated flux of 
IRAS 23365+3604 is less than 3a, the peak intensity at the 
velocity resolution of 23 km/s is more than 5mK, which is 
about 4a as shown at the top of Fig 2. We consider that CS 
5-4 emission is detected in IRAS 23365+3604 and the flux 
derived from Gaussian fitting was used. 

To avoid the system uncertainty between different tele- 
scopes, we have not included the CS 5-4 data detected 
with other telescopes for our study. The fit results for 14 
sources with CS 5-4 detected sources is: Log(LiR)(LQ) = 



0.938(±0.1) x Log(L CS )(K km s^pc -1 ) + 4.156(±0.834), 
with correlation coefficient of 0.939. NGC 2146 and NGC 
1068 have the largest deviations from the fit as shown in Fig- 
ure 3. NGC 1068 is known to harbor active galactic nucleus 
with non-negligible contribution to the total 24 /im emission. 
So the infrared luminosity within the CS 5-4 beam using the 
aperture correction of the 24^im MIPS image could be over- 
estimated, and thus NGC 1068 is likely closer to the best 
fit. Our CS 5-4 spectrum has similar line sha pe and veloc- 
ity ra nge to those obtained with IRAM 30m (|Martm et al.l 
l2009t ), but with total flux about twice that from IRAM 30M. 
Since the beam size of IRAM 30m, which is about 10" for ob- 
serving CS 5-4, is much smaller than the CS emission region, 
which is about 15" fr om high resolution o bservation of CS 
2-1 with IRAM PdBI (|Tacconi et al.lll997f ). So our observa- 
tion with HHSMT can pick up more extended emission than 
that with IRAM 30m. On the other hand, a s a nearby LIRG 
with only a possible low luminosity AGN (|lnui et al.l 120051 ; 
llsraej 120091 ). the infrared luminosity of NGC 2146 within 
CS 5-4 HHSMT beam should not be over estimated as that 
for NGC 1068. The large deviation should not be from the 
offset problem of beam switching mode, because the 2' off- 
set used in the observation is much far from the disk region 
with molecular gas based on CO mapping observation ^Israeli 
120091) . Because of limited velocity coverage and broad line 
of CS 5-4 (see Figure 2) with the full line width of about 
700 km/s, which is con stant with the result of HCN 1-0 
jGao fe Solomonll2004al ), the over subtracted baseline of CS 
5-4 may cause under estimation of the line flux, which can 
cause the large deviation. Since it is hard to do re-calibration 
of CS 5-4 and/or infrared emission in NGC 2146 and NGC 
1068, we will leave these two sources in this sample as solid 
detection and do fitting with original CS and infrared lumi- 
nosity. We also tried to do the same fitting but without NGC 
2146 and NGC 1068, which gave the result: Log(L IR )(L Q ) = 
1.037(±0.06) x Log(L CS ){K km s~ 1 pc- 1 ) + 3.29(±0.50), 
with correlation coefficient of 0.981. 

The relationship between CS 5-4 and infrared luminosi- 
ties presented in this work is consistent with close associ- 
ations between star formation and dense gas which is also 
suggested by previous studies of correlations between SFR 
and HCN 1-0 (ISolomon et al.ll 19921 : [Gao fc Solomonll2004bh . 
With energy from upper level to the grand level E up ~ 35K, 
which is much higher than the temperature (10 — 20K) of 
quiescent molecular cloud without star formation, and a lit- 
tle bit higher effective critical density than that of HCN 3-2 
which is about 10 times that for HCN 1-0 (|Evanslll999l ). CS 
5-4 is a good tracer of active star forming gas. Dense gas 
mass in galaxies estimated from individual types of dense 
gas tracers, such as HCN 1-0 and CS 5-4, may have system- 
atic errors due to different physical conditions (temperature, 
density, etc.) excitation mechanisms (infrared pumping in- 
stead of collisional excitation) in different galaxies. . The 
contribution of nuclear activity to the total infrared lumi- 
nosity can cause the over estimation of star formation rate, 
which can also introduce scatter to the correlation. Even 
though the line ratios of different dense ga s tracers have been 
found to be variable in different galaxies (lAalto et al. I l2002l : 
IGreve et al] l2009h . the relation between dense gas mass, 
traced by both HCN 1-0 (iGao fc Solomon! l2004bh and CS 
5-4 in this work, and the star formation rate traced by in- 
frared luminosity has an unity slope. So we suggest that 
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dense gas, which is directly linked to star formation, is very 
important for understanding star formation in galaxies, and 
the relation between star formation and dense gas should 
be similar with the observations of different dense gas trac- 
ers, otherwise that dense gas tracer should be affected by 
other properties, such as shock or infrare d pumping. Recent 
theoretical work (|Krumholz et al]|2009l ) showed that both 
the total gas and the molecular gas star formation law can 
not even be fitted by a single power law. Theoretical predic- 
tions of dense gas star formation law and dense gas fraction 
in different type of galaxies are needed for comparing with 
observational results. 



4 SUMMARY AND PROSPECTS 

With the observations of CS 5-4 using the HHSMT 10m tele- 
scope toward 24 local infrared bright galaxies including 20 
LIRGs/ULIRGs and 4 nearby normal galaxies, we detected 
CS 5-4 emission in 12 LIRGs/ULIRGs and 2 nearby normal 
galaxies. As the first extragalactic CS 5-4 survey with a sam- 
ple larger than 10 galaxies, we found a correlation between 
CS 5-4 luminosity and infrared luminosity with an almost 
unity slope. 

Since only 14 galaxies have CS 5-4 detection in our 
sample, it is necessary to enlarge the sample for CS 5-4 
observation. The existing (sub-)millimeter telescopes, such 
as HHSMT, CSO, JCMT, IRAM 30m, APEX, and also 
the LMT can be used for such observations. The nearby 
star forming galaxies can be observed with the smaller 
telescopes (HHSMT, CSO, JCMT, and APEX), while the 
ULIRGs/LIRGs should be observed with IRAM 30m and 
the LMT which have better sensitivity. We have started 
a project to detect CS 5-4 in more ULIRGs/LIRGs with 
IRAM 30m. Other CS transitions, such as CS 1-0, 2-1, 3-2, 
and also 7-6 can also be considered to determine dense gas 
properties (volume density, temperature, etc.) for galaxies 
with strong CS emission based on one transition observa- 
tion. CS lines in galaxies should also be a good choice to 
study star formation in galaxies for ALMA early science. 
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CS 5-4 in galaxies 

Table 1. CS 5-4 survey source list and galaxy properties 



Source name R.A. DEC D Rcdshift Log(L j R ) 

(,12000) (J2000) Mpc Lq 



NGC 6946 


20: 


34 


52 


3 


+60:09:14 


3.65 





000133 


10 


16 


NGC 891 


02 


22 


33 


4 


+42:20:57 


1. 13 





001761 


10 


27 


NGC 2146 


06 


18 


37 


6 


+ 78:21:19 


12.6 





002979 


11 


07 


NGC 7331 


22 


37 


04 


1 


+34:24:56 


14.5 





002722 


10 


58 


NGC 1068 


(12 


42 


10 


7 


-00:00:48 


16.1 





003793 


11 


27 


NGC 2276 


07: 


27 


1 1 


3 


+85:45:16 


33.8 





008059 


10 


81 


NGC 4194 


12 


14 


Oil 


5 


+ 54:31:37 


35.5 





008342 


11 


06 


CGCG 049-057 


15: 


13 


13 


1 


+07:13:32 


55.5 





012999 


11 


27 


NGC 992 


02 


37 


2. r . 


5 


+ 21:06:03 


55.7 





013813 


11 


02 


IRAS 17578-0400 


18: 


00 


3,1 


9 


-04:00:53 


60.0 





014043 


11 


35 


NGC 7771 


23 


51 


24 


9 


+ 20:06:43 


60.9 





014267 


11 


34 


MCG+ 12-02-001 


00 


54 


03 


6 


+ 73:05:12 


67.1 





015698 


11 


44 


NGC 1614 


1 


33 


59 


9 


-08:34:44 


68.1 





015938 


11 


6 


NGC 7469 


23 


03 


15 


6 


+08:52:26 


69.8 





016317 


11 


59 


NGC 828 


02 


10 


Oil 


6 


+39:11:25 


76.8 





017926 


11 


31 


Arp 220 


15 


34 


57 


1 


+23:30:11 


77.6 





018126 


12 


21 


NGC 2623 


08: 


38 


2 1 


1 


+25:45:17 


79.3 





018509 


11 


54 


VV 114 


01 


07 


47 


2 


-17:30:25 


86.1 





020067 


11 


65 


NGC 6240 


16: 


52 


58 


9 


+02:24:03 


106 





024480 


11 


85 


NGC 6090 


16 


11 


40 


7 


+ 52:27:24 


126.6 





029304 


11 


51 


Mrk 231 


12 


56 


1 1 


2 


+ 56:52:26 


184.0 





042170 


12 


51 


IRAS 17208-0014 


17: 


23 


21 


9 


-00:17:00 


187.0 





042810 


12 


39 


VII zw31 


05 


08 


15 


3 


+ 79:36:46 


236.1 





053670 


11 


94 


IRAS 23365+3604 


23 


39 


01 


3 


+36:21:09 


286.0 





064480 


12 


13 



a. Total Infrared Luminosity is from RBGs (Sanders 2003), and 
corrected with the distances. 

b. The sources from NGC 6946 to NGC 2276 have been corrected 
from total infrared luminosity from this table to the emission 
within the region of CS 5-4 beam, using MIPS 24 /an image, to 
compare with CS 5-4 luminosity. 



Table 2. Results of CS 5-4 line observation with HHSMT 



Source name 


I CS5-4 


Area(CS5-4) 


AV CS5-4 


On source time 




K km/s 


K km/s 


km/a 


Minutes 


NGC 6946 


1.299 ± 0.097 


1.40±0.19 


480 ±70 


45 


NGC 891 


3.692 ± 0.327 


3.64 + 0.50 


377 ±54 


20 


NGC 2146 


0.611 ± 0.060 


0.60+0.07 


315 ±54 


110 


NGC 7331 


< 0.626 ± 0.266 






15 


NGC 1068 


0.883 ± 0.180 


0.67 + 0.14 


130 ±29 


60 


NGC 2276 


< 0.72 ± 0.361 






15 


NGC 4194 


0.976 ± 0.298 






45 


CGCG 049-057 


1.100 ± 0.304 


1.11 + 0.43 


304+105 


45 


NGC 992 


< 0.138 ± 0.159 






30 


IRAS 17578-0400 


0.908 ± 0.133 


0.90 + 0.16 


283±50 


45 


NGC 7771 


0.419 ± 0.124 


0.46+0.16 


296±88 


60 


MCG + 12-02-001 


0.642 ± 0.060 


0.69+0.08 


368±41 


70 


NGC 1614 


< 0.088 ± 0.044 






140 


NGC 7469 


< 0.244 ± 0.167 






40 


NGC 828 


< 0.202 ± 0.124 






20 


Arp 220 


1.631 ± 0.128 


1.70 + 0.22 


467+68 


60 


NGC 2623 


< 0.147 ± 0.083 






100 


VV 114 


< -0.441 ± 0.249 






15 


NGC 6240 


0.993 ± 0.118 


1.06 + 0.18 


289±58 


80 


NGC 6090 


< 0.110 ± 0.160 






100 


Mrk 231 


< -0.390 ± 0.117 






45 


IRAS 17208-0014 


0.981 ± 0.330 


0.86 + 0.32 


141+65 


20 


VII zw31 


0.393 ± 0.109 


0.46+0.12 


341 + 81 


140 


IRAS 23365 + 3604 


< 0.309 ± 0.177 


0.32 + 0.09 


64+25 


00 



a. less— 4 is the flux integrated from the spectrum with the line 
width of CO TO from literature. 

b. Area(CS5-4) is the flux using gaussian fitting for sources with 
solid detection except for NGC 4194, which has the line shape 
more complicated than one gaussian component. The derived line 
width was listed as AVcss— 4- 
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Figure 3. The 14 galaxies with CS 5-4 detection are plotted as 
'o', while the other 10 sources are plotted with left arrow and 'o' 
using the 3<r of the integrated flux to estimate the upper limit of 
CS 5-4 luminosity. 



Figure 1. CS 5-4 spectra of Arp 220 (top) and NGC 1068 (bot- 
tom) observed with HHSMT. 
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Figure 2. CS 5-4 spectra of IRAS23365+3604 (top) and NGC 
2146 (bottom) observed with HHSMT. 



